Magnetization transport in spin ladders and next-nearest-neighbor chains 
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We study magnetization transport at high temperatures in several spin ladder systems as well as 
in next-nearest-neighbor coupled spin chains. In the integrable ladder considered we analytically 
show that the transport is ballistic in sectors with nonzero average magnetization, while numerical 
simulations of a nonequilibrium stationary setting indicate an anomalous transport in the zero 
magnetization (half-filling) sector. For other systems, isotropic Heisenberg ladder and spin chains, 
showing eigenlevel repulsion typical of quantum chaotic systems, numerical simulations indicate 
diffusive transport. 



I. INTRODUCTION 



Understanding transport in quantum and classical sys- 
tems from first principles has a long history. Perhaps the 
simplest question one can ask is what is the nature of 
transport in a given system, for instance, is it ballistic, 
in which case localized disturbances spread with time to 
a region whose size grows linearly with time, or, is it 
diffusive, in which case a region will grow as a square 
root of time. While the situation is clear in systems of 
non-interacting particles - in the absence of external scat- 
tering effects non-interacting systems are ballistic - for 
interacting systems (also called strongly correlated) such 
question proves to be very difficult to answer, even in the 
simplest conceivable models. 

A paradigmatic example is, for instance, a one- 
dimensional Heisenberg model^'^. Although being orig- 
inally devised as a toy model explaining high tempera- 
ture of a ferromagnetic transition in metals, nowadays 
its anisotropic XXZ version serves as one of the simplest 
strongly interacting quantum systems. Despite being 
solvable by the Bethe ansatz^ its nonequilibrium physics, 
in particular magnetization transport, is being debated 
for many years. The main obstacle one has to face trying 
to assess transport is a lack of efficient out-of-equilibrium 
tools, while on the other hand evaluating linear-response 
formalism using the Bethe-ansatz solution seems to be 
too difficult, except in the simplest case of zero temper- 
ature^. One exception are integrable systems for which 
the so-called Mazur's inequality^'^ can be used to bound 
a current-current correlation function. Mazur's inequal- 
ity states that if any of an infinite number of constants of 
motion that an integrable system possesses has a nonzero 
overlap with the current, the correlation function will not 
decay, thereby proving ballistic nature of transport. In 
the Heisenberg model this is for instance the case with en- 
ergy transport^ or magnetization transport away from 
zero-magnetization sector®. In the zero- magnetization 
(half-filling) sector though all known local constants of 
motion have, due to the spin-fiip symmetry, zero over- 
lap with the magnetization current^ and no conclusion 
can be reached based on the Mazurs's inequality. Re- 
cently, a new quasi-local constant of motion has been 
found enabling one to prove ballistic transport in the 



gapless phase |A| < 1^. For |A| > 1 no rigorous re- 
sults are known and one has to resort to numerical cal- 
culations making situation less clear, often even contro- 
versial. For A > 1 and high temperatures magnetiza- 
tion transport seems to be difFusivei2r2S, although more 
involved picture sometimes emergesi^iSS. Also at tem- 
peratures below the ground state gap theoretical argu- 
ments suggests diffusion^i. Experimentally the most im- 
portant, but at the same time also the most difficult for 
numerical studies, is the isotropic point. At isotropic 
point A = 1 some numerical investigations suggest an 
anomalous transporti^^— , see also Ref. 12^ for indica- 
tions of an anomalous behavior at the isotropic point, 
while other— ii^i2^"— indicate ballistic transport or are 
inconclusive. 

As the above example of an integrable Id XXZ spin 
chain shows, transport behavior can be rather complex. 
At first sight appealing conjecture that, due to constants 
of motion, integrable systems display ballistic transport 
seems to be incorrect. Namely, as mentioned, an easy- 
axis XXZ Heisenberg chain shows strong evidence for dif- 
fusion. That solvability does not necessarily imply bal- 
listic transport is indicated also by the existence of a 
solvable diffusive, albeit dissipative. Id model"^^. Non- 
integrable, in general quantum chaotic"^^ systems, are 
on the other hand expected to display diffusion. How- 
ever, as recently demonstrated^, there are exceptions to 
this rule. It has been rigorously shown that in a spe- 
cial class of an XX-type spin ladders (that class, for in- 
stance, includes the Hubbard chain) ballistic subspaces 
exist. Though probably exceptional, these counterexam- 
ples show that there are no general rules connecting, e.g., 
integrability or chaoticity, with the transport behavior. 
In light of this it is important to gather information on 
a transport type in different chaotic and integrable sys- 
tems. 

In the present work we shall study magnetization 
transport at high temperature in a number of spin ladder 
systems, including integrable and non-integrable ones. 
Transport, of either magnetization or heat, has been 
studied beforeiSiii^'SiiH^i^ for some of the systems we 
consider, for other the presented results are new. Note 
that next-nearest-neighbor coupling terms, frequently 
used to break integrability in spin chains, in fact result 
in ladder systems. Therefore, some of the ladder systems 



studied can also be considered as spin chains with next- 
nearest-neighbor (n.n.n.) couphng. We should also men- 
tion that spinfuU particle chain models like, e.g., the Id 
Hubbard model^ '^^i^^'^^i'^^' '— , can be, via Jordan- Wigner 
transformation, rewritten as spin ladder models. Spin 
ladder systems are not just of theoretical interest but 
are realized in a number of compounds, for a review see 
Refs. 46 andli^. 

Apart from one integrable ladder, we shall exclusively 
focus on systems with strong chaos and at transport 
properties at high temperatures and half-filling. We shall 
demonstrate that the magnetization transport is diffusive 
in such systems. In the integrable ladder we on the other 
hand find anomalous transport at half-filling and analyt- 
ically prove ballistic transport in sectors with non-zero 
magnetization. 



II. METHODS 

There are different ways to numerically assess quantum 
transport. One is via linear response theory by evaluat- 
ing equilibrium time-dependent current autocorrelation 
function. In numerical calculations one is always lim- 
ited to finite-size systems causing two effects: for finite 
L and time t the correlation function C{t) might not yet 
converge to its thermodynamic limit value and, going 
with i — >■ CO the correlation function will not decay to 
zero, even in a diffusive system, but will rather have a 
finite-size fluctuations. Therefore, due to these finite-L 
and finite-i effects great care is needed to correctly evalu- 
ate the limits limt_^oo limi^oo. Another way of studying 
transport is to directly simulate nonequilibrium states. 
There are two possibilities, one can study transient dy- 
namics of initial nonequilibrium states like, e.g., spread- 
ing of localized packets and calculating how fast their 
width increases with time, or, one can go to a stationary 
setting in which a constant driving is applied to the sys- 
tem. The latter approach has the advantage that there 
are no finite-time effects, only finite-size, as one by defi- 
nition studies a nonequilibrium stationary state reached 
after an infinite time. In present study we shall use a 
nonequilibrium stationary setting. 



A. Lindblad master equation 

Nonequilibrium situation will be induced by a bound- 
ary coupling to magnetization reservoirs. These can, with 
certain probabilities, fiip the boundary spin either up or 
down. If these probabilities are different at two ladder 
ends the driving will cause a nonequilibrium state. Spin 
fiips at the boundary are described in an effective way 
with the so-called Lindblad operators Lk, while the den- 
sity matrix describing the ladder evolves according to the 
Lindblad master equation, 

dp{t)/dt = \[p(t), H] + C^'^ipit)) = Cip{t)), (1) 



C'-ipit)) 






Lkp{tlLl] + [Lk,pml]- 



Provided the dissipative part C^^^ is nonzero one will typi- 
cally have a single stationary state poo, being the solution 
of £(poo) = 0, to which an arbitrary initial state p(0) con- 
verges after long time, poo = limt^oo p{t)- In a nonequi- 
librium setting such a state poo is called a nonequilibrium 
steady state (NESS). 

Before specifying in detail Lindblad operators used, let 
us comment on the applicability of the Lindblad equa- 
tion within the context of quantum transport. Lindblad 
equation can be derived from microscopic equations of 
motion of the system plus reservoirs under certain, from 
the condensed-matter perspective, rather restrictive con- 
ditions of a weak coupling and a fast decaying environ- 
mental correlationsi^i. While these conditions are usually 
well satisfied in, e.g., quantum optical systems, where 
the environment is fast, this is not so in condensed mat- 
ter. Environmental degrees there (electrons in the leads, 
phonons,...) are not necessarily fast compared to the 
timescale of the system of interest. As a consequence, 
evolution equation for the system will not be local in 
time, like the Lindblad equation ^, but will in general 
be non-local with a nontrivial integral kernel accounting 
for memory effects. While memory-effects can play a role 
in a transient finite-time behavior they are not expected 
to be important in the long time limit of nonequilibrium 
stationary states considered here. In certain situations 
one can even show exactly that memory effects, i.e., non- 
Markovian effects, play no role for NESS^. 

Lindblad operators modeling the reservoirs will differ 
depending on whether we want to study ladders or next- 
nearest-neighbor chains (that can be viewed as ladders 
with a diagonal inter- rung coupling, see Fig. [SJ. For 
ladders both spins in the first and in the last rung are 
coupled to the reservoir. Lindblad operators that we use 
are 



ii,2 - Vr(i T m) f^i , i3,4 = vr(i±70'^L: 

^5,6 = Vr(lTM)Ti±, L7^s^VW±V)ri, 



where a^ and r^ are Pauli matrices on the 1st and the 
2nd ladder leg, respectively, and cr^ = {a^ ± i(T^)/2, 
T^ = {t^ ± iT^)/2 (in Eq.([2]) upper and lower signs in ± 
go together). L is the number of rungs. For n.n.n coupled 
chains only the left-most and the right-most spins are 
coupled to reservoirs. Lindblad operators are 



Li. 



rViTM'^r, L 



3,4 



r^l±pa^. 



(3) 



For chains L is the chain length. The coupling F in both 
cases plays no essential role and we fix it to F = 1. The 
most important parameter is the driving p,. For zero 
driving p = the NESS is a trivial poo ^ 1, i.e., a equi- 
librium state at infinite temperature. For small /i the 
stationary state is a true nonequilibrium state, though 
it still stays close to an identity density matrix meaning 
that we are studying nonequilibrium systems at an infi- 
nite (high) temperature. Driving is also symmetric with 



respect to the left/right end and the NESS obtained has 
always zero expectation value of the total magnetization 
(i.e., we are at half- filling). Note that the choice of actual 
Lindblad operators is quite arbitrary. The choice used in 
the present work is perhaps the simplest and has been 
used in a number of our previous studies, see also, e.g. 
Ref. [50- Provided the boundary effects are small, which 
is the case at high temperature^^, and for non-ballistic 
systems, the bulk properties should be independent of 
the details of driving. 

Ladder and chain systems that shall be considered all 
conserve the total magnetization in the z direction. The 
corresponding unitary symmetry is f7 = exp (— ia X], '^'X)i 
with UHW — H . Because a dissipative Lindblad term 
C^^^ (|2l3p is also invariant under such [/, nonequilibrium 
steady states considered in the present work are all in- 
dependent of an optional homogeneous magnetic field in 
the z direction added to H . That is, if poo is the NESS 
for C with _ff, then the same p^o is the NESS also for 
C with H' = H + B J2j ^] ■ This is a general con- 
sequence of the symmetry of the master equation. To 
see this let us denote by 1^ a general unitary symme- 
try, and by C the corresponding conserved quantity. Let 
V be an exact symmetry of the Liouvillian ©, that is 
VC{p)V'^ = C{VpV'^). First, one notes that a unique 
NESS is invariant under V, VpooV^ = poo, see e.g., 
Ref. i55|- This means that in the eigenbasis of a corre- 
sponding conserved quantity C the matrix poa is block- 
diagonal, while the matrix C is diagonal with identical 
elements within each diagonal block. C and poo therefore 
commute and, if p^o satisfies i[poo,H] -K £^"'(/9oo) = 0, it 
also satisfies i[poo,H + C] +£^^^{pod) — 0, i.e., poo is also 
the NESS state for H' = H + C. 



B. Numerical method 

Because we want to study system's properties in a sta- 
tionary state we have to obtain poo- There are two pos- 
sibilities: one can either solve the stationary equation 
'C(poo) = 0, or, one can integrate the Lindblad equation 
^, obtaining p{t) and from it the NESS in the limit 
t -^ oo. We use the latter method by first writing p{t) 
in a matrix product form with matrices A^. of fixed di- 
mension M, describing a site k and an element s of a 
local operator basis. Ladders as well as n.n.n. chains are 
treated as a ladder system with an arbitrary coupling 
between two nearest-neighbor rungs. One rung is con- 
sidered as a single site k, so that the dimension of the 
operator basis at one site is 4^ (i.e., the number of dif- 
ferent values of the index s in matrices A|). The total 
number of complex parameters describing a state p{t) of a 
ladder with L rungs is therefore IQLIvP. Choosing large 
enough M any state p{t) can be written in such a matrix 
product operator form. Time evolution is then evaluated 
using the time-dependent density renormalization group 
method^ by writing a short-time propagator e'~' as a 
series of single and two-site transformations. The method 



we use is an adaptation^^ for dissipative systems in which 
the optimality of a matrix product decomposition is pre- 
served by reorthogonalizations, for details see Ref. [Sj. 
Evaluating two-site transformations exactly would lead 
to an exponentially increasing (in time) matrix dimen- 
sion M . Numerically this can not be handled and one 
truncates dimension after each transformation to a fixed 
size M . This truncation is the main source of errors in 
the numerical method. How large should M be depends 
on the amount of entanglement that a state p{t) has in 
the operator space. For instance, the equilibrium state at 
an infinite temperature is a product state (a product of 
identities at each site) and can be represented by matri- 
ces of size A/ = 1. For small /z, where p{t) is still close to 
1, we therefore expect that one can do with a reasonably 
small M . This is a reason why simulations of NESS at 
high temperatures require the smallest M . In our simu- 
lations we used matrix sizes of up-to M = 150 and ladder 
lengths of L < 100, requiring several weeks of simulation 
time for a single run on a 10-core CPU computer. 



C. Assessing transport 

Once the NESS is obtained - after time t that typically 
scales as an inverse power of L - expectation values of any 
local operator can be evaluated. Our main focus is on the 
magnetization current and on the magnetization profile 
along the ladder/chain. Fixing driving /i, typically at— 
p = 0.2, we study how magnetization current j scales 
with the system length L. If the scaling is j ^ 1/i, the 
system is diffusive and obeys a phenomenological trans- 
port law 



J = ~D^z, 



(4) 



where Vz is a magnetization gradient and Z? is a size- 
independent transport coefficient (diffusion constant). 
Other extreme situation would be when j is indepen- 
dent of i, signaling ballistic transport. Transport that 
is intermediate between ballistic and diffusive is called 
anomalous and is signaled by a scaling j ~ 1/L" with 
< a < 1. The nature of transport can be also inferred 
from magnetization profile. For small magnetization a 
profile is linear for diffusive systems while it is non-linear 
in the case of an anomalous transport where D can be 
considered to be length-dependent^^. 

All above statements about the scaling should be in- 
vestigated in the thermodynamic limit L -^ oo. Even 
though we limit ourselves, besides one integrable case, 
to fully chaotic systems, where the convergence with L 
is expected to be the fastest, it turns out that in some 
cases sizes L ~ 100, though rather large for a quasi-exact 
simulation of a strongly interacting quantum system, still 
might not be large enough to reach thermodynamic limit. 
Going to significantly larges sizes is at present not pos- 
sible due to rapidly increasing simulation times. Simula- 
tion time increases with L because the convergence time 
to Poo increases, but even more significantly, the required 



matrix size Af also increases because the observables, like 
the current, decrease with L and so a larger M is required 
to obtain the same relative accuracy in j . 



D. Checking for quantum chaos 

Despite exceptions^-, one in general expects that for 
non-integrable strongly interacting quantum system, in 
other words for systems displaying characteristic features 
of quantum chaos, transport is diffusive. In order to con- 
vince ourselves that the systems we study are indeed not 
being close to integrability, we have checked their chaotic- 
ity by calculating the spacing distribution of nearest en- 
ergy levels. The so-called level spacing distribution p[s) 
is a standard criterion of quantum chaos in Hamiltonian 
systems^. In chaotic systems there are no selection rules 
and the eigenenergies will "repel" each other, leading to a 
deficit of small spacings between two consecutive eigenen- 
ergies. In quantum chaotic systems with time-reversal 
symmetry the expected theoretical level spacing distribu- 
tion is well described by the so-called Wigner's surmise 
for the orthogonal ensemble. 



p{s) = s— exp (— s^7r/4). 



In integrable system there are selection rules due to con- 
stants of motion, resulting in an exponential form of 
p{s) = exp(— s). One should bear in mind that in or- 
der to see chaotic level statistics ([5]) explicit symmetries 
of a system have to be taken into account. Spacing has 
to be calculated within a single symmetry class. The 
symmetries of the systems studied are described in Ap- 
pendix |X1 



III. RESULTS FOR LADDERS 

There has been a number of works studying magneti- 
zation transpor t ^°i^^'^^i'^^ as well as heat transport^!"— 
in spin ladders. For studies of transport in Id Hubbard 
model, that can be equivalently rewritten as a spin lad- 
der, see Refs.i,[ni[ia[2ll,|3i-[4l The prevailing conclu- 
sion is that in non-integrable ladders at high temperature 
transport is diffusive, being in-line with a general rule of 
expecting diffusion in non-integrable systems. Numeri- 
cal studies have been mostly limited to systems with less 
than 1/ = 20 rungs; in present work we shall study signif- 
icantly larger systems. For a review on ladder systems, 
includin g re ferences to an extensive experimental work, 
see Refs-liel andlJTl 

We shall name different ladder systems according to 
the type of the coupling between nearest-neighbor sites. 



XX- type is a coupling of the form cr^cr^_ 



i+l 



'i+l 



XXZ- 



type is a coupling of the form a'fafj^i -|-crf tTf^_;^ -I- Aaf af,^! , 
while an isotropic Heisenberg coupling is equal to the 
XXZ coupling with the anisotropy A — 1, i.e., an XXX 
coupling. 



Ladder systems that shall be studied are depicted in 
Fig. [T^ and Fig. [T}d. The same methodology that we use 
has been used before to study the so-called XX ladder. 
Fig. [It, with an XX-type coupling along the legs and an 
XXZ-type in the rungs. A special case of such an XX lad- 
der is the Id Hubbard model obtained if the coupling in 
rungs is a'^rf . As shown in Ref . |44| the Id Hubbard model 
is diffusive under symmetric driving at infinite tempera- 
ture. This diffusive transport is not changed in the pres- 
ence of an additional XX-type coupling in rungs^^. It 
is instructive to rewrite a tight-binding system of free 
fermions on a ladder in spin language. Namely, for 
a tight-binding model we know that is is ballistic be- 
cause it is equivalent to a system of free fermions. Us- 
ing Jordan- Wigner transformation it can be written as 
a ladder shown in Fig. [TJi with a 4-site coupling inter- 
changeably connecting neighbors in the upper/lower leg 
being of the form (c^cr^+i +CTfcCrfc+i)T^T^^;^ (written here 
for the upper leg). We have numerically checked (data 
not shown) that such a ladder coupling indeed results in 
ballistic magnetization transport. Note that the r^r^,,.]^ 
term in the above coupling is absolutely crucial for the 
ballistic transport to appear; without it one would have 
an ordinary XX-type ladder displaying diffusive trans- 
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FIG. 1. Schematic representation of different spin ladders: 
a) integrable ladder, Eq. ^, b) isotropic Heisenberg ladder, 
Eq. ([6]), c) XX-ladder, d) free fermions on a ladder. Straight 
line denotes an XX-type coupling, spring a ZZ-type coupling, 
double line is an isotropic Heisenberg coupling while a straight 
line with two springs in d) is a coupling involving 4 sites (see 
text). Full points mark sites that are coupled to a reservoir, 
Eq. ©. 
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FIG. 2. Level spacing distribution for the isotropic Heisenberg 
ladder (|6]). U — 1, L — 8 and four symmetry sectors with 
Z = and zero total spin size are used (1026 spacings in 
total), see Appendix fXl 



A. Isotropic Heisenberg ladder 



Isotropic Heisenberg ladder is described by 



L-l 



^-E 



CTi ■ CTi+l + Ti ■ Ti+l 



^E' 



(6) 



Isotropic Heisenberg ladder, and in particular its ver- 
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chaotic as indicated by good agreement of level spacing 
distribution with the Wigner 's surmise ^ demonstrated 
in Fig. m 

Regarding magnetization transport, in Ref. 13 it has 
been found that initial localized packets spread out dif- 
fusively. In our stationary nonequilibrium setting we use 
a symmetric driving of Eq. © so that in the NESS 
Poo magnetization flows only along both legs, and is 
equal in the two legs, while there is no current in the 
rungs. The current operator is defined via a continu- 
ity equation for local magnetization ct| + ^k^ resulting 
in jl°^ = i[al + T^,hk,k+i], where /ifc,fc+i is the local 
hamiltonian density. For the model in Eq.([6]) we obtain 
JT*" — Jk + Jl ' with the current operator in the upper leg 
j^ — 2{a^(j^_^_^ — cr^cr^^j^), and a similar expression for 
jl in the lower leg. Due to a symmetric driving of both 
legs ([2]) in NESS both currents are the same and, due to 
continuity, independent of the site k. We shall therefore 
simply study the current in one of the legs and denote 
J — (Jk) ~ (Jk)^ with averages being expectation values 
in the NESS poo ■ In Fig- [3] we show current and magne- 
tization profile Zk = tr {a^poo) in the NESS. At /x = 0.1 
the results would be almost indistinguishable from the 
presented ones for /x = 0.2 in Fig. [3^. While magneti- 
zation profile is linear (with only very small deviations 
at few edge sites), suggesting diffusion, the scaling of the 
current shows small deviation from a diffusive ^ 1/ L. At 
U = 1.0 the scaling is j - l/i"-^, while at [/ = 0.5 it is 
j ^ 1/L°-^^. Note that at [/ = one would have two un- 
coupled isotropic Heisenberg chains for which an anoma- 
lous j ^ 1/i"-^ scaling has been observedi^i^. From the 
finite-size data presented it is difficult to make a defi- 
nite conclusion whether magnetization transport in the 
Heisenberg ladder is diffusive or not in thermodynamic 
limit. Considering a rather linear magnetization profiles 
we deem it plausible that small deviations observed are 
due to finite-size effect and the transport would become 
diffusive in thermodynamic limit. 
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FIG. 3. a) Scaled magnetization current on L, b) magnetiza- 
tion profile in one of the chains {U = 1,L — 64). All for the 
isotropic Heisenberg ladder ^ and ^ = 0.2; in a) data for 
U — 1.0 and U = 0.5 is shown. 

sion with different coupling strengths along rungs and 
legs (U ^ 1), is experimentally very much relevan t ^^'^^ . 
The system ([5]) has a nonzero spin gap^. It is quantum 



B. Integrable ladder 



The Hamiltonian is 



L-l 



•x,+i)+4t/^^,.r,. (7) 



It is a Heisenberg ladder with an additional four-spin in- 
teraction^. At U = the model is called a spin-orbital 
model^ and can be obtained as a large-f/ limit of a two- 
orbital Hubbard model at quarter fiUing'^^. Spin-orbital 
model can be, up-to an irrelevant constant, written as 
Hsu{4} = H{U = 0) = Y.J P],j+i, where P,j+i is a per- 
mutation operator on two rungs, Pjj_|_i|a,/3) = |/3,a), 
and \a), |/3) are two arbitrary rung states. Alternatively, 
it can be expressed in terms of generators G"'^ = W){M 
of the SU(4) group, i?su(4) ^ Efe Ej.,a G'fe ^^^+1- Inter- 
action is therefore SU(4) invariant and the spin-orbital 



model can be considered to be a generalization of the 
isotropic Heisenberg chain (that has an SU(2) symme- 
try) and is sometimes called a SU(4) Heisenberg model. 
It is Bethe ansatz solvable in one dimension by the gen- 
eral method^^ for systems with permutation interaction. 
One-rung operators Ci — J2j '^] +"^1' ^"2 = J2j '^j'''j a-^d 
C3 = J2j ^^'''f + '''J'^J + '''i'^l ^''^ conserved quantities 
for any fJ, while aX U — Q also all three components of 
^ ■ (Tj and X^iT? ^'''^ conserved. Nonzero rung interac- 
tion U ([7]) , being equal to C3 , plays the role of a chemi- 
cal potential, preserving integrability of the system^! for 
any U . The model is gapless^ for U <1 and gapped for 
U >\. 



1. Ballistic transport in nonzero-magnetization sectors 

Let us for a moment consider i?su(4) obtained for 
U — Q. Because the Hamiltonian is a sum of nearest- 
neighbor transpositions one can easily construct invari- 
ant subspaces that will display ballistic transport. Tak- 
ing the singlet \S) = (|01) - |10))/V^ and triplet states 
IT) = (|01) + |10))/y2,|O> = |00),|/) = 111) for a 
rung basis, and for instance the initial state of the lad- 
der \S . . . SIIIS ■ ■ -S), we can see that i?su(4) acting on 
such a state will cause the left-most and the right-most 
Is to spread ballistically to the left and right, respec- 
tively, causing two ballistic fronts. Although such con- 
struction is similar to the one in Ref. l33l two situations 
are fundamentally different. i?su(4) is integrable, and, 
as we shall show, energy current is a constant of motion 
causing ballistic transport in sectors with nonzero mag- 
netization, whereas the XX ladder discussed in Ref. [Sj 
is quantum chaotic with a more complicated dynamics 
than just transpositions (for instance, there the energy 
current is not a constant of motion). To see why the 
model ([7]) is ballistic away from zero magnetization (half- 
filling) sector let us first define currents. The magnetiza- 
tion current operator is independent of U and is j|,°* = 
j'^{l+Tk-Tk+i)+Jl{l+crk-(Tk+i), where j^'^ are the same 
chain currents as for the isotropic Heisenberg ladder. Lo- 
cal energy current, defined by jf = i[hk~i,k,hk,k+i], is a 
sum of an [/-independent term and a term proportional 
to U, jf = j|(f/ = 0) + f/ • jf{U ^ 0). jf{U = 0) 
is simply the energy current of the iJsu(4) model and is 
(written for fc = 2) 



jf{U = 0) = [{atj^, + al^l^ + al3l2) (8) 

with jjj^ = 2(cr^(T^-crJcr^) and /ij^ = TyTk, and (cr or) 
meaning all preceding terms with a and r matrices being 
interchanged. One can show that, taking periodic bound- 
ary conditions, the total energy current of the 7Jgu(4) 
model, J^ = J2kjfi^ = 0)j is an exact constant of mo- 
tion, [jf, H] = 0, regardless of U (note that J^ =J2k ^k 
however, is not). Because Jq has in addition a nonzero 
overlap with the magnetization current it can be used 



to bound the spin Drude weight away from zero. Let 
us denote by J^ = J2k^k°^ ^^^ total magnetization cur- 
rent. Thermodynamic overlaps needed, e.g., {JqJ^), are 
relatively simple to evaluate at infinite temperature but 
finite chemical potential, where the grand-canonical state 
is p = Yik Pk^ with pk ~ exp {—4><jX) being an equilibrium 
state of one spin. Identifying z = tr (p^cr^) as the aver- 
age equilibrium magnetization, or, equivalently, as the 
filling fraction / = (z + l)/2, all averages are polynomial 
functions of z. Denoting by Da, the magnetization Drude 
weight, Mazur's inequality^ can be used to obtain 



Ds > f i^, 



K = -- 



1 Jn^J' 



E tS\2 



K = T 



L {Jfjf) 
0^(l-z^)(l + z2)^ 
15 + 13z2 + 10z4 + 2z6 



(9) 



holding at high temperature. Observe that the bound is 
independent of C/, simply because Jq and J^ are, even 
thou the Hamiltonian in Eq. ([7]) does depend on U . Away 
from maximal polarization, 2: 7^ ±1, and nonzero mag- 
netization, z ^ 0, the value of K is nonzero, proving 
nonzero spin Drude weight in the integrable spin ladder 
([7]) at high temperature, and, as a consequence, ballistic 
magnetization transport. In present work we shall nu- 
merically study the case z = (/ = 1/2) where there 
remains a possibility to have non-ballistic transport. 



2. Numerical results for zero magnetization sector 

In numerical simulations we shall use the critical U — 
1. Note however, that, due to the fact that the interaction 
term U is equal to the constant of motion C3 of the closed 
system, magnetization transport for our symmetric driv- 
ing is almost independent of the value of U . If the sym- 
metry V = exp (— ia^j, <Tj; • Tfc) corresponding to the 
conserved quantity C3 = X]fe ""fe ' "nc, would be an exact 
symmetry of the LiouviUian ©, VC{p)V'^ = C[VpV'^), 
then the NESS state poo would be exactly independent 
of U . In our case the symmetry V preserves the unitary 
part, VHV'^ = H , but is not an exact symmetry of the 
dissipative part ^. Therefore, in an open system V is 
only an approximate symmetry; it is violated at bound- 
aries. Still, we find^ that the magnetization transport is 
almost independent of U . This also shows that the size 
of the ground state gap by itself does not play any role 
on transport at high temperatures. 

In Fig. m we show the scaling^ of magnetization cur- 
rent of one chain species j = {j'^{l+Tk-Tk+i)) with L and 
one instance of the magnetization profile. The current 
scales as j ^ 1/_L°-^^, indicating anomalous transport. 
Correspondingly, magnetization profile along the ladder 
is not linear but rather displays larger gradients towards 
the ends. Similar profiles have been observed^^'^^ in 
the isotropic Heisenberg model, also showing anomalous 
transport j ^ 1/L" with a = 1/2. Note that both, 
the isotropic Heisenberg model and the integrable lad- 
der ([7]), are special due to their SU(2) and SU(4) sym- 
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L-1 



(b) 




FIG. 4. a) Scaled magnetization current and b) magnetiza- 
tion profile in the upper leg {L — 64) for the integrable ladder 
system, Eq. ((T]). In b) the dotted line is the best-fitting gra- 
dient used in the scaling of current in frame a). U — 1, and 
driving /i = 0.2. 



metry, respectively. On a speculative note, considering 
that a = 1/2 for the SU(2) model, and a = 2/3 for 
the SU(4) one, general rule would be that the expo- 
nent of anomalous transport is a = N/{1 + N) for a 
model H ^ J^k Pk,k+i with an SU(2N) invariance. This 
would suggest diffusive (a = 1) transport in classical 
limit TV — >■ oo. For a recent study of transport in classical 
Heisenberg model see Ref.^. 



IV. NEXT-NEAREST-NEIGHBOR CHAINS 




FIG. 5. Different spin chains with a next-nearest-neighbor 
coupling: a) XX chain with a ZZ n.n.n. coupling, b) XX 
chain with an XX n.n.n. coupling, c) isotropic Heisenberg 
chain with a ZZ n.n.n. coupling (a straight line is an XX- 
type coupling, a spring a ZZ-type coupling while a double 
line is an isotropic Heisenberg coupling). Full points mark 
the sites that are coupled to a reservoir Q. 




FIG. 6. Level spacing distribution for spin chains, a) Isotropic 
Heisenberg with a ZZ n.n.n. coupling, b) XX chain with a XX 
n.n.n. coupling, c) and d) XX chain with a ZZ n.n.n. coupling. 
All data is for L = 8 and a sector with Z = (averaging over 
4 subsectors; in total around 11000 levels for each system), 
see Appendix!^ Full curve is the Wigner's surmize ([SJ. 



It is believed that integrability-breaking perturbations, 
provided they are large enough, will in general induce dif- 
fusive transpori*iai^i20i25^^-CT^ gome studies?5-6§ ob- 
served indications of ballistic transport at low temper- 
atures. In present work we reconsider the question of 
magnetization transport in spin chains with integrability- 
breaking next-nearest-neighbor coupling in the regime of 
strong integrability-breaking (quantum chaos) . Different 
spin chains studied are shown in Fig. [5j Note that by 
numbering ladder sites in a zig-zag manner, the next- 
nearest-neighbor coupling of a chain is in the ladder for- 
mulation given by the coupling terms in both legs, while 
the nearest-neighbor coupling of the chain is the ladder 
coupling in rungs and the diagonal inter-rung coupling. 



A. XX chain 

First, we are going to study the XX chain with a ZZ 
next-nearest-neighbor coupling. 



L-l 



L-2 



-I-2- 



(10) 



Magnetization current in the NESS is a standard j = 
(2(cr^cr^^^ — cr^o-fc^i)). Integrability breaking perturba- 
tion of strength [/ = 0.5 and U — \XS shall be used, for 
which the system is quantum chaotic. In Fig. |5] we can 
see nice agreement of level spacing distribution with the 



8 



Wigner's surmise. 

We shall also study the XX chain with a XX type n.n.n. 
coupling, 



L-l 



L-2 



iJ = 



= 1 



{afa^^,+4a\^,) + uTafa- 



=1 



1+2 



'i+2- 



(11) 

Ax, U — 0.5 studied the model is again quantum chaotic, 
see Fig. [6l Magnetization current operator gets in this 
case an additional next-nearest-neighbor term, and is 

J, = 2{alal^, - alal^,) + 2U[alal^^ - al^l^,). (12) 

As one can see in Fig. [7] magnetization transport is in all 
cases diffusive, indicated by the scaling j ~ 1/i, as well 
as by linear magnetization profiles. 



(a) 



^ 0.1 



0.01 




L-1 



(b) 




140 



FIG. 7. a) Current scaling for the XX chain with a ZZ 
(squares and circles) or XX (triangles) type n.n.n. coupling. 
Scaling is in all cases diffusive. Straight lines are 16.6/1/, 
9.0/1/ and 3.8/1/ (top to bottom), b) Magnetization profile is 
linear (L = 128). jx = 0.2. 



B. Isotropic Heisenberg chain 

We also studied the Isotropic Heisenberg chain with a 
ZZ n.n.n. coupling. 



L-l 



L-2 



iJ = ^ ^, • cr,+i + [/ ^ o^cr, 



i+2 J 



(13) 



and we shall use U — l.Q for which the system is quantum 
chaotic, Fig.ini Magnetization current is j = (2((t^(T^^-^ — 



CT^o-fc+i)). As shown in Fig. [8] the current scales as j ^ 
1/L^-^, while the profiles show slight deviations from a 
linear function close to chain ends. Note that, as is most 
often the case, if the current scales faster than ~ 1/L, 
i.e., if the system goes towards an insulating regime, a 
local gradient in profiles is smaller close to the system 
edge (Fig. ^, on the other hand, if the scaling is slower 
than ~ 1/L, i.e., if the system goes towards ballistic, 
the gradient is larger (e.g.. Fig. H]). In the present case 
deviations are small and it is difficult to asses if it is just 
a finite size effect and the system becomes diffusive in 
thermodynamic limit. 
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FIG. 8. The isotropic Heisenberg chain with a ZZ n.n.n. cou- 
pling of strength (7=1.0, /i = 0.2. a) scaling of current with 
system size, b) magnetization profile for L = 100. 



V. CONCLUSION 

We have studied magnetization transport in a linear re- 
sponse regime at high temperature by numerically calcu- 
lating nonequilibrium stationary states of Lindblad mas- 
ter equation. For the isotropic Heisenberg ladder, being 
quantum chaotic, we find close-to diffusive behavior, with 
the differences being possibly due to finite-size effects. In 
the XX spin chain with strong next-nearest-neighbor in- 
teraction transport is always found to be diffusive. The 
Isotropic Heisenberg chain with an integrability-breaking 
next-nearest-neighbor interaction is also very close to dif- 
fusive. We also found that the integrable ladder, which 
at t/ = has a SU(4) symmetry, shows anomalous mag- 
netization transport at half-filling, while away from half- 



filling, using Mazur's inequality, we prove that the trans- 
port is ballistic. 
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Appendix A: Symmetries 

We study magnetization transport, i.e., transport of 
the z-component of spin. All systems considered there- 
fore conserve the total magnetization in the z-direction. 
For ladders this is the operator Z = X],=i '^\ + '''fc- The 
corresponding symmetry transformation is a rotation 
[/z — n, exp(— iaapexp (— iarj). Under U Pauli ma- 
trices transform as U a'^XJ'^ — a^,, Ua^W = cos(2a)cr^ + 



sin (2a)cr^, and UalU^ = — sin (2a)(T| + cos {2a)al, and 
similarly for r^. 

There are two geometrical symmetries of the underly- 
ing ladder lattice. One is a parity P^ in the x-direction, 
obtained by mapping of sites /c — > L -f 1 — fc, while the 
other is a parity Py in the y-direction, obtained by the 
swapping of the two legs, cr^ ^^ ^k ■ 

In addition, there is a spin-flip symmetry given by 
the transformation U = J| '^l^f^ ^■^■' ^ rotation 
exp {—YKa'^/2) around the x-axis. It changes the sign of 
cr^'^ while it preserves ct^. It commutes with the rotation 
Z/z around z only in the sector with zero total magneti- 
zation Z = {). 

Symmetries of the isotropic Heisenberg ladder de- 
scribed by H in Eq. ([6]) are both parities P^ and Py, 
spin-flip and total magnetization Z . In addition, square 
of the total spin (^ . <7j -\- tj)^ is also a constant of mo- 
tion. 

For chains with a n.n.n. coupling, Eqs. (|10lllll3p . the 
symmetries are spin-flip, total magnetization Z and a 
product of parities PxPy 
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